Extended molecules and geometric scattering resonances in optical lattices 
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We develop a theory describing neutral atoms scattering at low energies in an optical lattice. We 
show that for a repulsive interaction, as the microscopic scattering length increases, the eftective 
scattering amplitude approaches a limiting value which depends only on the lattice parameters. In 
the case of attractive interaction a geometric resonance occurs before reaching this limit. Close to 
the resonance, the effective interaction becomes repulsive and supports a weakly bound state, which 
can extend over several lattice sites. 



Recent advances in cooling and trapping of cold atomic 
gases (see 0] for a review) have led to the experimen- 
tal realization of optically trapped lattices gases 0, 1^ ■ 
For a dilute gas and at low energies all physical prop- 
erties of such gases can be expressed in terms of the 
two-body scattering amplitude It is well estab- 

lished that confinement of a gas in one or more spatial 
dimensions can strongly modify the coUisional proper- 
ties of atoms, and can even induce geometric resonances 
0, 0, 0, • In an optical lattice related phenomenona 
should occur when, (i) the microscopic free-space scat- 
tering length becomes comparable to the size of a sin- 
gle lattice well, a situation readily obtained by combin- 
ing Feshbach or shape resonances with optical trapping 
techniques [a llfl llll 113 ; ^^'^ (ii) the tunneling be- 
tween the neighboring lattice sites is sufficiently small. 
Below we give a theoretical analysis of the modification 
of neutral atoms scattering at low energies in an optical 
lattice. This serves as the basis of writing down an ef- 
fective (long wavelength) field theory of a dilute optical 
lattice gas with an effective atomic mass and effective 
atomic scattering length, modified and tunable by the 
lattice parameters. 

Let us consider neutral atoms confined in a 3D optical 
lattice and interacting through a short range potential. 
Our goal is the derivation of the low energy limit of the 
effective scattering amplitude, and in particular the ef- 
fective scattering length in a tight binding model 
(TBM) as a function of the free space scattering length 
a and lattice parameters, adapting techniques developed 
in 0] . We will show that the wavefunction ipp of the 
pair of atoms averaged over many lattice sites satisfies 
the effective Schrdinger equation for the relative motion, 



^^aeff<^(r-)^|r| ) Vp(r) = €pi^p{r), (1) 



where tp is the relative energy, and mg is the band mass. 
The main results are: (i) For repulsive interaction, when 
the microscopic scattering length a increases, the effec- 
tive scattering amplitude remains finite and continuously 
approaches a universal limit, which only depends on the 
optical lattice parameters, (ii) In the case of an attractive 
interaction a geometric resonance in the effective scatter- 
ing amplitude, similar to the resonances found in 01^, 
occurs before the universal limit is reached. In addition. 



after crossing the resonance the effective interaction be- 
comes repulsive and weakly bound states result close to 
the resonance. This predicts the existence of large lattice 
induced molecules which can extend over several lattice 
sites, which - by analogy with the Wannier-Mott exci- 
tons in semiconductors [lj| - will be referred to below as 
Wannier-Mott molecules. The parameters effective mass 
TOs and scattering length OeS appearing in are the rel- 
evant parameters for the long wavelength effective field 
theory describing the dilute lattice gas. 

We turn now to a discussion of two-body scattering 
in an optical lattice. The particles interact via a short 
range interatomic potential U{r) depending only on the 
relative coordinate r. The range of this potential is 
assumed to be small compared to the lattice spacing and 
to the size of the ground state wavefunction in an in- 
dividual well. In addition, the lattice spacing is assumed 
to be small relative to the de Broglie wavelength of the 
incoming particle. The particles move in an 3D optical 
lattice created by counterpropagating laser beams. The 
resulting periodic potential is characterized by a depth Vq 
and lattice spacing d, which are assumed to be equal for 
all spatial direction. Close to these minima, the potential 
can expanded as 
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where SR and Sr are the deviation of the center-of- 
mass and the relative coordinates R = {xi + X2)/2 and 
r — X2 — xi from their potential minimum positions. 
Here, xi and X2 are the positions of the particles and 
^ — m/2 and M = 2m are the reduced and the total 
mass, respectively. Moreover, the curvature within a sin- 
gle lattice site is given by = Vo/m(P. From here on 
we will use units such that h = m = 1. 

In such a periodic potential, the two-body wavefunc- 
tion 'I'(xi,X2) is characterized by certain quasimomenta 
Q and q of the center of mass and the relative motion, 
as well as by band indices S and s. Using the fact that 
the lattice potential is approximately separable (compare 
Eq. (0) ) , we can expand the two-body wavefunction in a 
set of basis functions 0, , 
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where W and w are Wannier functions. Within the TBM 
the functions W and w coincide with the stationary states 
of the Hamiltonian (jJJ. 

The wavefunctions of the relative motion (Bloch 
waves) are characterized by the quasimomentum q, the 
band index s and the energy esq = eg — ts cos{qd), where 
ts = y/DsUj/n, and Ds 1 is the WKB tunnehng expo- 
nent between the neighboring wells, and are the ener- 
gies of states s in isolated wells Q . For sufficiently small 
momenta the dispersion relation can be approximated as 
Esij ~ Cfi + q'^/'m-sj where = 2/tsCp is the effective 
mass. In the TBM, Zq < d, the tunneling coefficient is 
small and hence 3> 1. The asymptotic form of the 
relative part of the two-body wavefunction far from the 
region around r = where the interaction potential acts, 
can be written as 
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where fss' is the s— wave scattering amplitude, and g^' = 
^Jms'{e.sp — Cs'o)- This definition ensures that in the con- 
tinuum case / — ATTa/rris, where a is the free space two- 
body scattering length and rus — m. 

Since the interparticle interaction U{r) is very short- 
range, we can use free (Bloch) solutions for the relative 
motion everywhere apart from the immediate vicinity of 
r = 0. To calculate the effective scatteringamplitude 
we closely follow the method suggested in 0, la]. The 
most general solution of the free Schroedinger equation 
for 7' 7^ for a given energy e, which also contains the 
incoming wave is 

*(r) =^ii;,(r-r„)e^(P-'^")+^G(r,0), (5) 



The coefficient A can be found by considering the short 
distance assymptotics of the solution Indeed, at 

small distances r ^ Iq the wavefunction of the relative 
motion should match the solution of the two-body scat- 
tering problem in 3D ,15] : 
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where i? is a constant. The scattering length a takes 
into account all the processes occurring at distances of 
order tq <C I n, d - In particular a may contain the effect 
of Feshbach [1^ or shape resonances , provided that 
the size of the molecular bound state is sufficiently small. 

We consider incoming particles at low energies so that 
the scattering occurs in the lowest Bloch band s — Q. 
Taking the limit both r, r' ^ in the Eq. © we find 
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As it is well known the Green function is singular when 



r and can be represented as [IE 
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where T is the regular part of the Green function. The 
singular part comes from the direct classical trajectory 
connecting the points r, r', which has a universal char- 
acter. The regular part consists of the contribution of 
recurrent trajectories and depends strongly on the de- 
tails of the confining potential. 

The only resonance in the denominator of Eq.Q oc- 
curs for the s = contribution. Wc should therefore 
study this term separately by writing G = Gq + G' , where 
the function Gq is not singular, and 



where A is an arbitrary constant. Here the first term 
describes the incoming wave, whereas the second term 
gives the scattered waves. The Green function is given 
by its usual expression 



s' ,g,n,n' ^ ^ 

(6) 

For sufficiently large r we can expand esq up to second 
order around q = and integrate over q 
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where V is the volume of the system. If e^'o > ^sp then 
the corresponding partial wave does not propagate to 
r oo, qs' is imaginary, and the channel s' is closed. 
Comparing Eqs. Q,® and we identify the scatter- 
ing amplitude 

fss'ip.q') = A{esp)Vws{Q)5ss', 



Go = E 



(11) 



The TBM wavefunctions can be represented as Ws = 
ZsTpsif) where ips is the wave function of the oscillator in 
a state s, and Zg is a normalization factor (Z^ — N~^, 
where N is the number of the lattice sites). Substituting 
the wavefunctions in Ea. (|ll|l we find for the imaginary 
part of Go : 



Im Goip) = -t'^d^\MO)\^P- 
An 



The real part is given by the principle value 



Re Go(p^ 0) = - 
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The remaining terms in Eq.® can be analyzed in 
the WKB approximation. In a spherically symmetric 
harmonic trap the spectrum of the excitations is = 



3 



^eff 








I* 



Figure 1: Dependence of the effective scattering length a^g 
on the absolute value of the microscopic scattering length \a\. 
The dashed line corresponds to repulsive interaction potential 
(a > 0). The solid line corresponds to the attractive potential 
(a < 0) with a geometric resonance at \a\ = h. 



Lo{nx + Uy + Uz + 3/2). The WKB wavefunctions in each 
direction look like 



_ u; cos($„(a:)) 



where $n(a; — > 0) ~ 7Tn/2 + qx with q = \J^(n + 1/2). 
Then 
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where the summation occurs over even values of n. Note 



that the term with rix = ny 



is absorbed in Gn 



and hence is excluded from the sum. This sum diverges 
when X = x' . Changing from summation to integration 
and separating out the singular contribution wc find: 



G'(r^ 0,0) = 
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The sum over p can be estimated to give a contribution 
~ to T. Within the TBM {Dq <C 1) this is much 
smaller than the contribution given by ReGo and can 
therefore be neglected: 
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Substituting Eq. in Eq. and comparing the re- 
sult with the boundary condition ^ we can compute 
the coefficients A and B to find the expression for the 
scattering amplitude: 
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where /3 = {d/lo)^, and the effective scattering length is 
defined as 
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with U = /o£'yV41n2. 

The dependence of the effective scattering length QcS 
on the microscopic scattering length a is depicted in 
Fig. 1. First, according to Eq. (|15l) in the limit of small 
a << Z,, we have acs/a — j3ms >> 1, i.e. the effective 
interaction in the lattice in comparison with free space 
is strongly enhanced by both the strong confinement of 
atoms in the lattice sites f3 = [d/lof » 1 and the 
large effective mass rris . In the case of repulsive interac- 
tion, a > 0, Oeff continuously grows as a increases and 
approaches the universal value = msPl* ~ d > Iq. 
In the case of attractive interaction, a < 0, the effective 
scattering length exhibits resonant behavior and diverges 
at \a\ — l^. This is the same type of geometric resonance 
suggested for quasi-lD and quasi-2D gases 0j 11] at 
low energies. At larger values of \a\ the effective scat- 
tering becomes repulsive and reaches the universal limit 

Ooo ■ 

The imaginary part in Eq. (|14|) is required by the opti- 
cal theorem and does not imply any inelastic process 
(apart from those possibly contained in Ima) . Depending 
on the sign a^g the poles of the scattering amplitude give 
rise either to bound or to virtual states in the effective 
potential as a result of the interplay of the interparticle 
interaction and the lattice potential. An interesting pos- 
sibility arises when the microscopic scattering length a 
is negative and its absolute value is slightly larger than 
the resonance value l^. In this case the effective interac- 
tion is repulsive and is characterized by a very large and 
positive scattering length 



\a\ — L 



Accordingly, the scattering amplitude H14(l has a pole at 
\p\ = P* corresponding to a bound state with a binding 
energy 
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The size of such a dimer is ^ ~ and must be very 
large: ^ » d so that the bound state of the two parti- 
cles can extend over many lattice sites, and thus should 
be referred to as "Wannier-Mott" molecules in analogy to 
the Wannier-Mott excitons in semiconductors jl3| , where 
bound states of electron and holes can very large com- 
pared with the lattice constant. 

Note that in the TBM l.^. < Iq and, therefore, the 
geometric resonance \a\ ~ can occur at realistic val- 
ues of the microscopic scattering length a. Because the 
rate of 3-body recombination scales as a^g 19| the spon- 
taneous formation of these molecules should be facili- 
tated in the limit that the effective scattering length be- 
comes large. However, at the same time, in this limit the 
binding energy is very small and a reasonable fraction 
of bound atoms can only be expected at temperatures 
T ^ |e* |. A possible way of detecting these extended and 
weakly bound Wannier-Mott molecules arises if they are 
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composed of fermionic atoms in different internal states. 
When the trapping potentials for the separate states are 
tilted in opposite directions, the unbound atoms experi- 
ence a force and will move in one or the other direction, 
depending on their internal state, whereas the molecules 
experience zero net force. If the potential gradient differ- 
ence exceeds a certain critical value related to the binding 
energy, the molecules dissociate and the remaining atoms 
are released. 

The limiting value of the 2-body scattering amplitude 
/oo(0) is of particular importance for condensed matter 
calculations 
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and determines the effective interaction strength U = 
/oo(0) for a dilute lattice gas. For example, in the case 
of a Bose-condensed lattice gas the chemical potential is 
given by Un, where n is the density of the condensed par- 
ticles 20j . Similarly, in a system of interacting fermions 
with an attractive interaction (U < 0) the critical tem- 
perature of the BCS transition is 



Tr 



i exp(7 — 2) 



ep exp(-27rVpF|f/|), 
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where ep = Pp/'^Tns is the Fermi energy. An increase of 
the absolute value of the scattering length leads to a res- 
onance in the effective interaction potential and causes 
the latter to change from an attractive to a repulsive in- 
teraction at \a\ = U. Note however, that close to the 
resonance the imaginary and energy dependent part of 



the scattering amplitude cannot be neglected and Eq. H16() 
loses its validity. A further increase of |a| leads to a re- 
pulsive interaction which destroys the BCS ground state 
and results in the formation of dimers. Since the lat- 
ter are bosons, at sufficiently small temperatures a BEC 
of dimers can be created. This possibility is related to 
the ongoing discussion of the BCS-to-BEC crossover in 
strongly interacting Fermi-gases with attractive interac- 
tion (221 . We note that the BEC of large molecules pro- 
vides an example of a superfluid ground state for a sys- 
tem of interacting fermions with effective repulsion. This 
possibility is related to the reported BEC of excitons in 
semiconductors [2§j. 

In summary, we have derived an analytical expression 
for the scattering amplitude as a function of the free space 
scattering length and the parameters of the confining op- 
tical potential for small energies. For a repulsive interac- 
tion, as the microscopic scattering length increases, the 
effective scattering amplitude grows and continuously ap- 
proaches a universal limit. In the case of attractive in- 
teraction, increasing the scattering length leads to a ge- 
ometric resonance occurring before the universal limit is 
reached. Close to the resonance, when the effective in- 
teraction becomes repulsive, it supports a weakly bound 
Wannier-Mott molecule which can extend over several 
lattice sites. 
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